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Abstract. - When non- magnetic impurities are introduced in a d-wave superconductor, both 
thermodynamic and spectral properties are strongly affected if the impurity potential is close to 
the strong resonance limit. In addition to the scalar impurity potential, the charge carriers are 
also spin-orbit coupled to the impurities. Here it is shown that (i) close to the unitarity limit 
for the impurity scattering, the spin-orbit contribution is of the same order of magnitude than 
the scalar scattering and cannot be neglected, (ii) the spin-orbit scattering is pair-breaking and 
(iii) induces a small id xy component to the off-diagonal part of the self-energy. 



In high-r c superconductors, disorder has important effects on both thermodynamic and 
spectral properties. The critical temperature T c and the superfluid density p s are lowered 
by non-magnetic impurity substitution ^, ^, ^ and disorder induced by irradiation ^| . 
Recent ARPES data show clearly how disorder leads to a redistribution of spectral intensity 
by adding new states at the Fermi level M. The basic elements of the current theory have 
been inspired by previous studies on heavy fermion superconductors and are given by the 
anisotropy of the order parameter and the strong resonance limit for the impurity potential 
[^, M. These elements, adjusted to describe condensates with a d-wave symmetry of the order 
parameter, are able to account for most of the features observed by experiments on high-T c 
d-wave superconducto rs [p| , [Tc| , Oj. However, discrepancies still exist, like the overestimation 
of the T c -suppression ]I2|. In order to correct this situation, and to provide a more realistic 
picture, several improvements of the theory have been proposed [H], Oj, and, recently, the 
effect of spatial variation of the order parameter has been taken into account jl2], [la, . 

In addition to the scalar impurity potential, the charge carriers are also spin-orbit coupled 
to the impurities. So far, this additional scattering channel has not been considered because 
the spin-orbit interaction is believed to provide, if any, only negligible effects (at least in 
the absence of a Zeeman magnetic field). This argument is based on the observation that 
the spin-orbit potential is of order v so ~ Ag v where v is the impurity potential and Ag is 
the shift of the g-factor |lS| . The value of Ag depends on the specific impurity, however 
its order of magnitude is roughly Ag ~ 0.1. From this estimate, it is expected therefore 
that the spin-orbit scattering rate l/r so ~ N v^ o , where N is the charge carriers density of 
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states, should be at most of order 1/t so ~ 10 _2 /Ti mp , therefore negligible with respect to 
the scalar impurity scattering rate l/ri mp JTg[ ] . Although such an estimate is correct in the 
Born approximation (weak scattering) nevertheless it underestimates the effect by orders of 
magnitude in the strong resonance limit, believed to be valid for high-T c superconductors. To 
illustrate such a substantial discrepancy between the Born and the unitarity limit, let us first 
consider the self-consistent t-matrix solution of the impurity problem for Ag — 0. In this case 
the electron (hole) propagator is G _1 (k, iu> n ) = i£b n — p^e{k) — piA(k), where A(k) = A cos(20) 
is the cZ-wave order parameter and <j) is the polar angle, e(k) is the electron dispersion for an 
half-filled band and p\, p% are Pauli matrices. The renormalized Matsubara frequency iuj n 
satisfies the following equation ]|, : 

iuj n = tuj n + — (i) 

cr - go(iVn) 

where go(iu> n ) = {i& n /[&n + A(k) 2 ] 1 / 2 } and (• ■ ■) is the average over the polar angle 4>. In 
equation (|l|), Y = rii/^No) and c = \/{ttNqv) where m is the impurity concentration. 
For c ^> 1, eq.(|l|) reduces to the Born limit while for c -C 1 it leads to the unitarity, or 
strong resonant, limit. Now, let us suppose that the spin-orbit impurity scattering leads to a 
renormalization contribution of the same form of eq. ([!]). Since v so /v ~ Ag, the renormalization 
induced by both the impurity and the spin-orbit scatterings can be estimated by: 



- go(iuJ n ) 2 (c/ Ag) 2 - g Q (iui n ) 2 



lUJ n =iaJ n +T — jj ■ ^ + ^ /n/A ^ 2 77— "T2- (2) 



For c > 1, iu> m ~ iu> n + T gQ(iuj n )c~ 2 (1 + Ag 2 ) and, as expected in the Born limit, the con- 
tribution of spin-orbit scattering is Ag 2 times smaller than that of scalar impurity scattering. 
On the other hand, for strong scattering, c/ Ag can be very small and eventually it vanishes 
in the unitarity limit c — > 0. As a result, in this limit the spin-orbit scattering leads to a 
renormalization of the same order of the impurity scattering, namely of order T / A. 

To provide more solid grounds to the above simple picture, it is necessary to treat the 
impurity and spin-orbit interactions on the same level by generalizing the usual t-matrix 
approach for the impurity scattering also to the spin-orbit contribution. To this end, let us 
start by considering the spin-orbit impurity potential. Two-dimensionality is often assumed 
in describing the main electronic excitations at least for some of the high-T c superconductors. 
In the present context, the reduced dimensionality has the following implication. If the 
charge carriers are confined to move in the x-y plane and the impurity potential is V(r) = 
v^2 i j i exp[ik • (r — Rj)], where denotes the impurity positions, the spin-orbit interaction 
assumes the following form: 

V so (r) = ^r, so ^e ^k ■( I ■- R ^k x p] • cr = % Ag v^^-^^^a^ (3) 

where p = — iVr is the momentum operator and the spin-orbit coupling has been parametrized 
by r] so v — Agv/kp where kp is the Fermi momentum. The effect of two-dimensionality is 
therefore to couple the electron spin only along the z direction. Hence, by choosing the z 
axis as the direction of spin quantization, the spin-orbit interaction (^|) does not mix the spin 
components. The generalized Green's function in the particle-hole spin space is: 

G(k, iu n )- x = G (k, iu n )- x - S(k, iw n ), (4) 

where Go(k, iw„) _1 = iuj n — P3e(k) — p2T2A(k) and the Pauli matrices pi and tj act on the 
particle-hole and spin subspaces, respectively pOI]. In the self-consistent t-matrix approxima- 



C. GRIMALDI: SPIN-ORBIT SCATTERING IN ETC. 



3 



tion the self-energy is S(k, iLu n ) — n.iT to t(k, k, iLu n ) where the t-matrix satisfies the following 
equation: 

T tot (k,k',^„) = u(k,k') +5^«(k > k")G(k" J iw n )r to t(k" > k',iw w ), (5) 

k" 

where u(k, k') = p 3 v + ir 3 Agv[k x k'] z is the impurity potential including the spin-orbit 
contribution. Because of the angular dependence of the spin-orbit interaction, it can be shown 
that the t-matrix (^) reduces to T tot (k, k'iw„) = T{iui n ) + T so (k, k', iuj n ) where 

T(iuj n ) = p 3 v + p 3 v ^ G(k, iuj n )T(iuj n ) (6) 
k 

is the usual momentum-independent contribution from non-magnetic impurities [g] and 

T so (k,k',^„) = ir 3 Agv[k x k.% + iAgvJ2$ x k"] z T 3 G(k", iu n )T so (k", k', iw n ), (7) 

k" 

is the t-matrix for the spin-orbit interaction. Before proceeding with the complete solution 
of eq.(^|), it is useful to analyze the lowest order contributions in Ag. By replacing G(k, iuj n ) 
with Go(k, iuj n ), the expansion of eq.(Q) up to the third order in Ag leads to a spin-orbit part 

(2) (3) 

of the self-energy of the form £ so (k, iu> n ) — Eso (k, iuj n ) + Y^ so (k, iui n ), where the first term 
is the usual Born contribution H^J (k, iuj n ) — niAg 2 v 2 J2k' 1^ x k'| 2 T3Go(k', iuj n )T 3 which 
renormalizes both the frequency and, contrary to the normal impurity scattering, the gap 
function. The term proportional to Ag 3 is instead: 

Si 3 )(k,^„) = -m^Agvf ^[kx k%[k' x k"] z [k" x k] z T 3 G (k', ^„)r 3 G (k", iu> n )r 3 

k',k» 

= — r A f/ - Sin ^ 2 (8) 
lW c 3 \ v /A(0') 2 +^/\v / A(0') 2 +^/ P U 

which contributes to the off-diagonal part of the self-energy. This term is proportional to 
i sin(2(/>) and has therefore a id xy symmetry which is a consequence of the angular dependence 
of the spin-orbit potential. Note that eq.(||) closely resembles the off-diagonal contribution 
found by Balatsky in the context of spin-orbit coupling of the condensate to magnetic scattering 



centres 21 . Although as a function of <ft, T, so (k, ioj n ) has maximum contribution where the 
gap function A cos(2</>) vanishes, it does not open any gap in the excitation spectrum. In 
fact, after analytical continuation iuj n — > lu, it can be shown that at the nodes of Acos(2</>), 
<fr = ±7r/4,±37r/4, the real part of the pole of the Green's function for k = kp satisfies 
lu = sio* (kp, lu) oc r(A<7/c) 3 (w/2A) log(2A/|o;|), so that lu — is the solution and no 
additional gap is opened. 

Equation (|J) suggests that the Green's function solution of the t-matrix problem is of the 
form: 

G(k, iLUn)' 1 = ioj n - p 3 e n - p 2 T 2 [A n ((f>) + ir 3 £i n ((p)], (9) 

where (D n , e„, A n (</>) and r2 n (</>) are frequency dependent quantities which must be calculated 
self-consistently. Note that in eq.(||), fl n (<p) is the off-diagonal contribution which reduces to 
eq.(^) at the lowest order in Ag. The t-matrix for the spin-orbit interaction (^) is solved by 
setting T so (k, k', iLu n ) — iAgv\k x t(k , iLu n )] z T 3 , where 

t(k',iw„) =k' + iAgv^2k"T 3 G(k",icu n )[k" x t(k' ,icu n )] z . (10) 
k» 
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The above equation is actually a system of two coupled equations for the components t x (k' , iuo n ) 
and t y (& ,iu} n ) which can be explicitly solved and, after some algebra, the spin-orbit contri- 
bution T so (k, k, iuj n ) to the total t-matrix becomes: 

T so (k, k,iu n ) = iAgvk x k y [Ayx(iuJn) ~ A~y(iu n )]T 3 

+ (Agv) 2 Y, [A^iiUnXkK^ + ^^nXkk,) 2 ] r 3 G(k',ia; n )r3,(ll) 
k' 

where A~y(iu) n ) and A~ x (iuj n ) are the inverse of the 4x4 matrices 



A xy (iu n ) = l~(Agv) 2 



A yx {iuj n ) = l-(Agv) 2 



^(k x ) 2 T 3 G(k, iuj n ) 
k 

^(fc,) 2 r 3 G(k, 
k 



^(fc,) 2 r 3 G(k, lUJ n ) 
. k 

lUJ n ) 

k 



(12) 



(13) 



Note that the first term in the right hand side of eq.(|ll|) is proportional to k x k y = (1/2) sin(2^>) 
and, in fact, it reduces to eq.(|J) at the lowest order in Ag. Moreover, this term would absent 
be for s-wave superconductors since in this case A X y(iui n ^ — Ay X (iuj n ^ . 

Plugging equation (j^) into ( |l2] ) and (|i"3| ) permits to invert A xy (iu) n ) and A yx (iuj n ), and the 
final equations for u> n , e n , A n ((f>) and f2 n (0) are obtained by demanding self-consistency for 
equations (^), (||), (^), (p4]-p^|). A consistent solution requires that A„(0) = A„ cos(20) and 
^n(</0 — f^ n sin(20), where 



A n 

S2j ; 



A + r 
-2r 



(c/Ag) 2 -(/ 2 -g 2 ) 
[{c/AgY - {gl + fl)f - Aflgl 
f n g n (c/Ag) 



fni 



where 



9n 



[(c/Ag)*-(g»+ft)]*-4ftg>' 
iu) n sm((j)) 2 \ / A n (4>) sm(4>) 2 



Jn 



'A B (0)a + a,(0) 2 + u>„ I \ V A„(0) 2 + £l n {4>) 2 + , 

while the equation for the renormalized frequency uj n is: 

, ^ g 0n , ^ (c/A 5 ) 2 + (/ 2 - 5 2 ) 



iui r , 



9dn 



{{c/AgY-igl+fDY-Aflgl 



fjn 



and 



90r, 



lA^y + n^y + coi 



(14) 
(15) 

(16) 

(17) 
(18) 



The set of equations (|T^)-(|T^) represents the main result of this paper and several features 
can be already outlined. First, although the equation for the renormalized frequency ( |l7| ) is 
more complex than the simple minded eq.(^|), the conjecture discussed in the introduction is 
confirmed, i. e., as long as c/ Ag <C 1 the renormalization due to spin-orbit interaction is of 
the same order of magnitude of the one induced by the impurity scattering. In addition, the 
spin-orbit scattering renormalizes also the gap function and induces the additional off-diagonal 
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self-energy Cl n . From eq.(|l5|) it can be seen that Q n oc (r/A)(A<?/c) 3 for Ag/c <C 1 and 
(l n oc (r/A)c/A<7 for c/ Ag <C 1 and is therefore much smaller than the spin-orbit parts of A„ 
and uj n which are of order (r/A)(Ag/c) 2 for Ag/c <§C 1 and T/A for cj Ag <C 1, respectively. 

Equations (|l4|)-(|l8|) must be completed with the equation for the order parameter A which, 
if the pairing interaction is V pa i r (k, k') = —V va ir cos(20) cos(2<//), reduces to: 



^T^^cos(20)Tr [p 2 r 2 G(k,^)] = AttT^ / hM^M ) , (19) 




n k 71 

where A = Vpair-No is the coupling constant. The critical temperature T c is obtained from 
cq.(|l9|) by setting A„ — > and f2„ — > in equations (|l4[ )-(|l9|). The resulting critical 
temperature satisfies a Abrikosov-Gorkov type of relation |22| : 

*(£)-G)-G^). 

where ^ is the di-gamma function and T c q is the critical temperature for the pure system 
(r = 0). The impurity and spin-orbit scattering parameters are 

r r - r 3(2c/Ag)2 + 1 (21) 

" c2 + l' so [(2c/A ff ) 2 + l] 2 - 1 j 

Superconductivity is destroyed for the critical value r* + T* Q = TrT c0 /2j (7 ~ 1.781). For 
Ag = 0, the c — > limit leads to the usual result for scalar impurities T* — ttT c0 /2j |ll| while, 
for Ag ^ 0, limc^oT* = kTcq/Aj, i. e., half the value of the impurity concentration in the 
unitarity limit without spin-orbit scattering. This result clearly shows that, in the unitarity 
limit, spin-orbit scattering has a strong pair-breaking effect of the same size of scalar impurity 
scattering. This feature is a consequence of the infinite resummation of the t-matrix which 
gives rise to a non-analytic point at c = 0, Ag = 0, and the limit Ag — > 0,c = differs from 
c — > 0,Ag = 0. The reduced critical temperature T c /T c q (solid lines) and the reduced zero 
temperature order parameter A/Aq (dashed lines) are shown in fig. 1 as a function of T/Aq for 
c = 0.1 and Ag — 0.0, 0.05, 0.1, 0.15 and 0.2 (from right to left). A is calculated numerically 
from the T — > limit of eq. ( p"9] ) and Ao is the gap for the pure case (in the weak coupling 
limit). As a function of impurity concentration, the suppression of both T c /T c q and A/ Aq is 
stronger for smaller values of c/ Ag. In the inset of fig. 1, T c /T c q and A/Ao are shown for the 
limiting cases c — > 0,Ag — (a) and Ag — > 0,c = (b). 

Another thermodynamic quantity of interest is the superfluid density p s which is calculated 
from 

Ps = 2nmN vlT V / A " Wc0S ^ )2 V (22) 

^\[A n (^ + a i (0) 2 +^]3/ 2 /' 

where vf and m are the electron Fermi velocity and mass, respectively p3]. In fig. 2 it is 
shown T c /T c q as a function of the zero temperature limit of p s /p s o, where p s o is the superfluid 
density for the pure system, for T/Aq = 0.1, Ag — 0.1 and different values of c. Note that, 
by lowering c, the curves move upwards and that this effect is visible also for the limiting 
cases c — > 0,Ag — (a) and Ag — > 0,c = (b) reported in the inset. This feature is given 
by the spin-orbit renormalization of A n , otherwise absent when Ag — 0. In Ref. it has 
been stressed that when T c /T c q is plotted against the zero temperature limit of p s /p s o the 
experimental data lay above the theoretical curve corresponding to the unitarity limit. As 
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inferred from both fig. 2 and the inset, the inclusion of spin-orbit scattering tends to cure this 
discrepancy. 

In addition to thermodynamic quantities, the spin-orbit interaction affects also the spectral 
properties of d-wave superconductors. In fact, nonzero values of Ag contribute to the amount of 
gapless excitations already provided by the resonant scattering with impurities. This feature 



can be investigated by performing the analytical continuation of eqs. (|l4)-(18) by setting 
iLu n — > lo and iuj n — > Co and then plugging the results into the quasiparticle density of states 
N(io) given below: 

N{UJ) S gn(£)Re ( _ ) . ( 2:! ) 



Nn 



/ Co 2 - A(<j>) 2 - n(4>) 2 , 



The resulting density of states is plotted in fig. 3 for T/Aq = 0.05, c = 0.05 (a) and for 
r/Ao = 0.1, c = 0.1 (b). In both cases, nonzero values of Ag provide an additional contribution 
to the gapless states lowering at the same time the intensity of the coherence peaks at w ~ A. 
In the insets of fig. 3 it is also shown how the zero energy density of states N(0)/No is affected 
by the spin-orbit coupling. 

In summary, it has been shown that, if the impurity potential is close to the unitarity 
limit, the spin-orbit coupling to the impurities is as important as the scalar impurity potential 
and its effects cannot be neglected. This result points toward a critical re-examination of 
the existing experimental data in addition to the recent theoretical developments based on the 
spatial variation of the order parameter. As a final remark, it should be noted that a promising 
experimental route for the estimation of the spin-orbit interaction in d-wave superconductors 
is provided by the Zeeman response under an external magnetic field [[l9| p4| . Preliminary 
results suggest that nonzero values of Ag drastically affect the Zeeman response of ei-wave 
superconductors p5|. 



The author would like to thank P. Fuldc and I. Vobornik for fruitful discussions. 
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Fig. 1. - Reduced critical temperature (solid lines) and order parameter (dashed lines) as 
function of r/A for c = 0.1 and Ag = 0.0, 0.05, 0.1, 0.15, and 0.2 (from right to left). 
Inset: the same quantities for the c — > 0,Ag = limit (a) and for the Ag — > 0,c = limit (b) 




Fig. 2. - Reduced critical temperature as a function of the reduced superfluid density for 
T = 0.1, Ag = 0.1 and different values of c. Inset: superfluid density in the c — > 0,A# = 
limit (a) and in the Ag — > 0,c = limit (b) 
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Fig. 3. - Quasiparticle density of states for T = 0.05, c = 0.05 (a) and T — 0.1, c = 0.1 (b) and 
different values of Ag. Insets: corresponding zero energy density of states as a function of Ag. 



